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ABSTRACT 


In this report a time-domain Green’s Function Method for unsteady supersonic poten- 
tial flow around complex aircraft configurations is presented. 

We focus here on the supersonic range wherein the linear potential flow assumption is 
valid. In this range the effects of the nonlinear terms in the unsteady supersonic compress- 
ible velocity potential equation are negligible and therefore these terms will be omitted in 
this report. 

The Green’s function method is employed in order to convert the potential-flow dif- 
ferential equation into an integral one. This integral equation is then discretized, in space 
through standard finite-element technique, and in time through finite-difference, to yield 
a linear algebraic system of equations relating the unknown potential to its prescribed 
co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex 
aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hy- 
perboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed 
to vary linearly within each panel. Consistent with the spatial linear (first-order) finite- 
element approximations, the potential and co-normalwash are assumed to vary linearly in 
time. 

The long range goal of our research is to develop a comprehensive theory for unsteady 
supersonic potential aerodynamics which is capable of yielding accurate results even in the 
low supersonic (i.e., high transonic) range. 
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LIST OF SYMBOLS 

speed of sound 

or 

contravariant base vector, see equation (25) 

«2 

contravariant base vector, see equation (25) 

$ 

(Mob - 1) 1/2 

E(P) 

domain function, see equation (15) 

Hp) 

Dirac delta function 

fijk. 

Kronecker delta 

Fik 

finite-element shape function 

G 

Green’s function 

H 

Heaviside function 

k 

reduced frequency, 

l 

reference length 

Moo 

free stream Mach number U^/a^, 

n 

unit normal to oy 

N 

unit normal to E# 

N e 

total number of elements 

N n 

total number of nodes 

u 

circular frequency 

P 

point having coordinates x,y,z 

P * 

control point, (x*,y*, 2 :*) 

r 

point having coordinates X, Y, Z 

p* 

control point, (X*,Y*,Z*) 

p-f- 

Hadamard finite part 

R 

P - P. 

R' 

hyperbolic radius, see equation (11) 

t 

time 

T 

nondimensional time, a^fH (l 

Uoo 

velocity of undisturbed flow 

x,y,z 

space coordinates 

X, Y, Z 

nondimensional Prandtl-Glauert coordinates 

a 

X = x/f3l,Y = y/l,Z = zfi 
sxirface of body in x,y,z space 

E 

surface of body in X , Y, Z space 

E, 

surface of element i in X, Y, Z space 

4> 

perturbation velocity potential 
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<fr nondimensional perturbation velocity 

potential, 0/C/ oo £ 

xjj co-normalwash in x,y,z space 

co-normalwash in X, Y, Z space 




Operators 

a 2 , a 2 , a 2 
^ 

Laplace operator in the physical space 

^ a 2 _ a 2 _ a 2 
ax 2 ay" 2 a/ 2 

ft 4 . n <1. 

at ^ u °°d x 

Supersonic dot product, see equation (5) 


i. INTRODUCTION 


In this report we demonstrate how the Green’s Function Method of Potential Aero- 
dynamics may be implemented in the time domain so as to enable it to handle unsteady 
supersonic flow around complex aircraft configurations. 

We focus here on the supersonic range wherein the linear potential flow assumption is 
valid. In this range the effects of the nonlinear terms in the unsteady supersonic compress- 
ible velocity potential equation are negligible, and therefore these terms will be omitted in 
this report. 

The Green’s function method (Ref. 1) is employed in order to convert the potential- 
flow differential equation into an integral one. This integral equation is then discretized, in 
space through standard finite-element technique (Refs. 2 and 3), and in time through finite- 
difference (Ref. 4.), to yield a linear algebraic system of equations relating the unknown 
potential to its prescribed co-normalwash on the surface of the aircraft. The arbitrary 
complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized 
into hyperboloidal (twisted qiiadrilateral) panels. The potential and co-normalwash are 
assumed to vary linearly within each panel. Consistent with the spatial linear (first- 
order) finite- element approximations, the potential and co-normalwash are assumed to 
vary linearly in time. 

A frequency-domain first-order Green’s function formulation for linear oscillatory su- 
personic flow has been developed recently (Ref. 5). Presented here is the corresponding 
theory formulated in the time-domain setting. 

The long range goal of our research is to develop a comprehensive theory for unsteady 
supersonic potential aerodynamics for the low supersonic (i.e., high transonic) range. 
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1.1 A Brief Description of the Green’s Function Method 

Before getting into the specifics of this report we begin with a brief description of 
the Green’s Function Method. This method applies to the equation of the perturbation 
velocity potential. The potential function $ at any point P* in the flow field is given by an 
integral of terms containing the value of the potential and its co-normal derivative on the 
surface, a, surrounding the body and its wake. An integral equation for the potential on the 
surface of the body is obtained by letting the point P t f approach a point on the surface. 
With this method, the wake is a natural by-product and is treated as a layer of doublets. 
It may be noted that the integral equation does not require that the boundary condition 
on the co-normalwash be satisfied, but rather makes use of the continuity of the potential 
as the control point approaches the surface a. The tangency boundary conditions are 
automatically satisfied by the type of representation obtained with the Green’s Function 
Method. 

In current applications, the surface of the aircraft is divided into small quadrilateral 
elements. Each element is replaced by a paraboloidal hyperboloid surface defined by the 
four corners of the element. In this process the continuity of the surface is maintained 
but discontinuities in the slopes are introduced. The aircraft wake, on the other hand, is 
divided into strips parallel to the streamlines. These wake strips originate from the trailing 
edge and extend to infinity downstream. It should be noted that integrals over these wake 
strips may be carried out in an analogous way to their subsonic counterpart (see Refs. 3 
and 6). 

In the 0 th order theory, the unknown $ (in the Prandtl-Glauert Space) is assumed 
to be constant within each element, while in the 1st order theory $ is taken in the form 
$ = $ 0 + -f- 77$ 2 + where £ and 77 are local element-wise surface coordinates, 

and the coefficients # 0 , . . . , $3 are chosen to interpolate the $ values at the four corners 
of the element. In either situation the integral equation is approximated by a system 
of algebraic equations. This system of algebraic equations is then solved by standard 
numerical methods. It has been found (see Ref. 7 ) that in the supersonic range at least a 
1-st order theory is required in order to yield a nonsingular set of algebraic equations due 
to a numerical rather than physical anomaly. The numerical implementation presented in 
this report employs l s ^ order panels. 


2. UNSTEADY SUPERSONIC FLOW 


Our point of departure is the linearized equation for the unsteady potential compress- 
ible aerodynamic flow 




n 2 

u oo 


d 2 (f> 
dt 2 


0 



where V„ is the Laplace operator in the physical space while 4 > is the perturbation potential. 
Choosing a frame of reference such that the undisturbed flow has velocity U ^ in the 


t The bar is used herein to indicate vector quantities. 



direction of the positive x-axis, the linearized total time derivative is given by 


d 0 


a 


dt at + u °° ax ■ 


( 2 ) 


Introducing the generalized Prandtl-Glauert transformation 

X = x/fil, Y = y/l , Z = z/L T = $ = i/U^l (3) 

where l is a characteristic length, M ^ — U^/a^ and /? = (A/Jb - l) 1 / 2 , Eq. (l) yields* 


£2$ 

VoV't’f P jj T - 2 + IMooqxWT ~ 0 

where o stands for the sxipersonic dot-product defined as 


(4) 


dob = ajbj — o,yby — (5) 

where a and b are two arbitrary vectors. Thus the operator Vo V stands for 

<9 2 d 2 c/ 2 

dX 2 dY 2 dZ 2 1 j 

Note that a o a is not necessarily positive. We define the ‘super-norm’ of a vector a as 

a' = ||a|| = |a o a| 1/(2 (7) 

and will use this notation in later sections of this report. 


2.1 Supersonic Integral Equation 

In order to obtain the Supersonic Green’s function integral equation we proceed as 
follows: With P and P t representing the sending and receiving points respectively, the 
Green’s function G for Equation (10) satisfies 

(9 2 G a^c 

v O V G + a\ ()T \ + 2M 00 = HP ~ P*i T — 1\) (8) 

G = 0 at oo 

one well known solution of which is given by (see Ref. 7) 


G 


H 

4t rW 


(<W + ^©-) 


where 


H(P, p .) = { 1 0 


if X - X t < 0 and R o R > 0 
otherwise 


* See Appendix C for derivation. 


( 9 ) 

( 10 ) 
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Here we define 


while 

with 


R = P — P* = (X - X t )i + (Y - Y t )j + (Z - Z t )k 
r' = PH - |(x - X *) 2 - (Y - y *) 2 - (z - z*) 2 l 1/2 

- A(T - T, + 0 J ) 


( 11 ) 


( 12 ) 


e ± = M 0O (X* - X) ± R' (13) 

Here {P\H(P, P*) = 1} defines the zone of influence, or Mach forecone, with vertex at P*. 


Next for a closed bounded surface E bounding a volume V, we define the domain 
function 

V 

V 


E (P* ) “ 1 lf — 

1 *' 10 if P* G 


(14) 


Note that for P* on E the function E(P t ) will measure the so-called supersonic solid angle 
of E at P* (see Ref. 7 for details). Hence E(P t ) satisfies the notation 


E(P,)=£Nov(^,y S + l 


(15) 


where N is the outward unit normal to E. Applying the Green’s Function Method, with 
the Green’s Function G given by Eq. (9), it can be shown (see Ref. 7 for details) that Eq. 
(8) may be integrated to yield 


4^p(p*)$(p*,r t ) - 0 ([fl/] e+ + [\T>]®“) 


-H 

IF 


dE + 




e 


d (H 


dE 


<9$ 

© + ae>+ 


©" 

dG~ 1 

[dT\ 

dN c 

dT 


dN c j 


; dE 


where the co-normalwash 'F is defined as 

f = V$o]V = d$/dN c 

with the conormal N c defined as 

N C = ~Njl + Nyj + N z k 

so that 


d 


-N r 


d 


d 


, 4 . TV 4- N 

dm ~ z dX^ y dY^ z dZ 


Furthermore 


0 ± = M^X - X t ) ± R 1 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
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It should be noted here that E includes the surfaces of the aircraft and its wake. 
Throughout this report, for simplicity, the condition of supersonic trailing edge will be 
assumed. Thus, in the absence of the wake, E represents the aircraft surface. However, 
the general method employed here is not limited to the supersonic trailing-edge condition. 
The treatment of the wake when trailing edges are subsonic is analogous to the time-domain 
subsonic formulations and hence the reader is referred to Ref. 4 for further details. 

3. NUMERICAL FORMULATION 

In this section, a space and time discretization procedure will be introduced in order 
to approximate the integral equation by a linear algebraic system of time-delay difference 
equations. Solving this time-dependent system in a step-by-step manner yields the desired 
time-dependent perturbation velocity potential solution on the aircraft surface. Once the 
velocity potential is known at a given time step, the pressure coefficient may be computed 
through Bernoulli’s Theorem. 

3.1 Finite Element Formulation 

Assuming that the surface E is divided into N e small finite elements E z , Equation (16) 
yields 

4«e(p.)*(p.,t.) = e{j (Jm B * + 

i - 1 ’ ' 

+ E// E .(W e+ +W e 4^(|)^ (2i) 

1=1 1 

^4 rr /rd$i e+ d0+ rd<&i e ~<90-\ f-H\ 

+ f-* Jh- U dT J Bn 7 + at ~dW) \~W ) dh 

i=i 1 

Each surface element E t is approximated by a hyperbolic paraboloid given in the form 

p = p c + Pi i + Piv + ? 3 eh (22) 

where P c ,P l ,P 2 and P 3 are obtained in terms of the locations of the four corner points 
as (See Fig. 1) 



Here, P + _, for instance, refers to the element corner for which £ — 4-1 and i] = —1. The 
other corners, P ++ , P_ + and P ._ are defined similarly. It may be noted that the surface 
defined, according to Eq. (22) is continuous since adjacent elements have in common the 
straight line connecting the two common corner points. 
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3.3 First Order Space Discretization 

In what follows we shall take the potential function &(£,??, r) over a surface element, 
say E ? , at a fixed time r, as 

$*(£,»?,*■) = [(i + OU - 0 $ +-( r ) + (i + 0(! + 0$++(0 

+ (1 - 0(1 + v)*-+(r) + (1 - 0(1 - 0 *— ( 0]/4 

( 28 ) 

Equation ( 28 ) expresses the values of at any point (£, rj, r) of Ej in terms of the values 
of $ at the four corner points-also at time r. 

More generally, consider the first-order global shape function with the following defi- 
nition: 

' (1 + £)(1 — 77) /4 if node k coincides with corner +— of element i 

(1 + £)(1 + 77) /4 if node k coincides with corner -+•+ of element i 

Fiki^ v) — (1 + 0(1 + 0/4 if node k coincides with corner — -f of element i 

(1 — 0(1“ 0/4 if node k coincides with corner of element i 

0 otherwise. 
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( 29 ) 


With Fm defined by Equation (29), Equation (28) may be rewritten as 

N n 

= '£F ik (t,r l )* k (T) (30) 

h~\ 

|i 

where N n is the total number of nodes on the surface £, and § k denotes $ at the k 111 
node. 

Similarly, the supersonic co-normalwash $(= N o V$) may be represented by (See 
Ref. 5) 

N n 

= '£F lk (d,r l )* k (T) (31) 

k ) 

The same first-order finite-element approximation, Equation (29), has been employed for 
if. Note that if # is approximated by the l' s * order finite-element expression while 
is represented by a O^ 1 order formula, ^*(^,r?,r) = \l>*(0,0,r) = const., a mixed type 
formulation would result. In subsequent portions of this report, $ and \& will be taken 1 st 
order. 

3.4 Numerical Approximation of the Integral Equation 

Due to the hyperbolic nature of the governing unsteady supersonic velocity potential 
equation, the Heaviside function appears in the resulting supersonic Green’s function inte- 
gral formulation, (eqs. (16) and (21)). The appearance of the Heaviside function H(P,P t ) 
under the integral restricts the integration to that portion of the panel that is within the 
Mach Forecone with vertex at P t . 

Note here that the potential, $, and the co-normalwash, if, under the integral are 
both confined to the Mach Forecone region. A Taylor expansion in both space and time 
may be written for # and $ so that they become explicitly dependent on the space and 
time variables. The spatial derivatives of the Taylor expansion for $(or\)/) may be com- 
puted through the first-order finite- element representation introduced in the preceding 
subsection. The temporal derivatives, on the other hand, may be approximated by a back- 
ward finite-difference in time. With $( 1 F) expressed as an explicit function of time, the 
constraint that $(^) varies within the Forecone may be imposed, through the Heaviside 
function, in a straightforward manner. 

We remark that the time variable includes a delay term which is a function of the 
spatial variables (see Eq. (13)). The kernels of the integrals are essentially of the type 
( R')~ n (where n equals 1 for source and n equals 3 for doublet). We adopt the tactic of 
treating R 1 , in addition to £ and rj, as independent variable. Thus, we express $ and if 
via Taylor expansion in these variables. Explicitly, let us rewrite the time delay expression 
as 

T - 0 ± = r(C,r/,T) ± R' (32) 


1 1 



where 


r{t,n,T)=T-M 0O \X<ti,ri)-X{t,ri)] (33) 

Thus, $ may be expanded in £,77 and R' around a point “o” which is the geometric center 


of the portion of the element that is in the Mach Forecone, i.e., 

± R') ~ $U 0 ,r, 0 ,T 0 ± R' 0 ) + ^ ^ 


(£ - Co) 


V "dr, dr dr. 


(r,-r , 0 )±^ | (£'-£') 


(34) 


/ d 2 $ d 2 $ dr dr d 2 $ dr dr\ & 

+ { d(dri + d(dr Wr) + dr/dr d£ + d^dlJr, ) J " ~ ~ 7,0 > 

+ sf‘f W *'-«) 2 

where the subscript o denotes evaluation at the geometric center of the portion of the 
element that is in the Mach Forecone while 



~~ = + x z n) 

(35) 

and 

= -'M.y, + X 3 () 

(36) 

with 




QP _ 

Xl = d c • * = (^1 + P 3V) • * 

(37) 


dp _ ^ _ 

= -g- ■ i = (P 2 + PaO ■ •' 

(38) 


d 2 P . . 

x ' d ~ d^d,, 1 ~ r * 1 

(39) 


Expansion (34) is 2 n<t order accurate. Consistent with the finite-element approximation 
for the space variables, the 2 nd order spatial derivative terms, i.e., d 2 <f>/d £ 2 and d 2 $/dr, 2 , 
and higher order spatial derivatives, vanish and hence it is sufficient to include up to 2 nd 
order terms only in Eq. (34). 


Next, recall in Section 3.3 that 


Nn 


**'(£, >7 ,r) = Eiq* (£,??) <S>*(r) 
k= 1 

where (£, rj, r) is given by Eq. (29) and r is a prescribed time. Hence we may compute 
the spatial derivatives of $ at the point o as 


fl$*(£o.» 7 o,r 0 ± R'o ) _ ^ dF ik {£o ,rio) ±R ,^ 




k 1 


(40) 
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(41) 


Making use of 
tion (21) becomes 


where & j k and ©y 
o and 


5n = E ail *‘( r » ± *»> 

a^'(i o ,ri 0 ,r o ±R' o ) ^a 2 F ik (( 0 ,n c )^ , „ M 

*) — " = L —gp%— ± R °) 

aH'{Urio,r„±R'o) _ SF, k (i 0 ^ a )a^ k (T 0 ±K) 

did t fa dr 

^ 2 ^(C„^or 0 ±7?') = ^ dFikjj g, do) d $>k{T 0 d : R' 0 ) 
dr) dr dr) dr 


(42) 

(43) 

(44) 


Equations (40) through (44) with P t coinciding with the node j, Equa- 


4t r£y$y(T) 

+ 5 y t |* t (r-e+) + * t (r-e; t )] 

+ Cj k [$ k [T - 0 + ) + * t (r - 0 " t )l (45) 

+ c tk \i k (T - etj + iqr - e^)] 

+ C ik [<i> k (T - e} t ) + i k (T - e>j k )] 

+ D, t {* k (T - e) t ) + i k (T - 6- k )} 

+ *it[*t(r-e+ t ) + i t (r-e; t )]} 

, are understood, for notational simplicity, to be evaluated at the point 


= E{E4*t(r-0+) + **(r -©,*)] 

• k=\ 1 

+ B jk \y t (T - e+ t ) + * t (r - ©jjj 


rr d F° d F° 

Bj t = E J] e \f& + gf (s - £o) + (n - -)„)+ 


B 


Jk 


d 2 F°, 

dWr)^° ~ ^o){d ~ 


-H 


R' 


d'Z 


+ 


= E/i i K s f[(s-s„) + ct| 

J(£ - Co)(?7 


dr 

<9t7 


3*1 dr 


di? fc dr| 

5?? di \ 


id - do) 

H 


R> 


dZ 


-i (if - i?;)(-A)< iS 


B, t = Y.Jj u (R' - R'o?{^t)^ 


(46a) 


(466) 

(46c) 



jpo f) / H \ 

+ ^ (e ” £ ") (r/ - r/ " )] ^U)^ (47a) 

<V= 

+ 

+ 

C,k = ^jfjR'-R'o) 2 j}f c ( §,)<*£ (47c) 

J= 1 1 

Af e /•/■ /_tf\ fiftl 

D ]k = ^ IL . ( ~J? T ) d]v T ( 48fl ) 

*“-1 1 

«,* = £ # E (*' - m 

t=l 1 

where F^, for instance, indicates F,*- evaluated at (£ 0 > *?<>)• 

Equation (45), with the coefficients defined by Equations (46) through (48), is the time- 
dependent aerodynamic operator relating the potential to its co-normalwash at N n nodal 
locations. These integral coefficients will be evaluated by a hybrid analytical-numerical 
scheme. More explicitly, each of these integrals will be integrated analytically in the in- 
direction followed by a Gaussian Quadrature numerical integration procedure in the r) 
direction. Discussion of these integrals will be presented in Section 4. 

In Equation (45) * the spatial variables have been discretized by a l s< order finite- 
element approximation. Following the procedure introduced in Ref. 4, discretization of 
the time-dependent functions, $ and will be discussed in the following. 

First, the time variable may be written as 


rr r)r r>T 

5 Ik ^ M |« ~ ^ + F ' k dv L<* _ ^ 


(R 1 - R, 


d t-H\ 
o) ON c \ R< y 


T — (n + a) AT where n is an integer and 0 < a < 1 (49) 

Next, a continuous function g(T) with T given by Equation (49) may be approximated 


g(T) = g[(n + a)AT] = (1 - a)g(nAT) + ag[(n + l)AF] 


H 



Setting T = nAT and using Equations (46) through (50), Equation (45) may be written 


as 


N n ___ , 

4tt Ej$j{n) = Y, - a1- k )V t (n - m+ ) + (1 


i)*k( 


k= 1 




+ B jk 


jk)^\^- <x jk)' i! k{n-m jk ) 

Jt - 1 )+° i ]k^k{n-rnJ k 

- m }k ) + ( J - - m 7*) 


n — m 


»} 


+ - m ik - 1) + «■,■***(« - m 




+ Bjk 

y { 
*jk 




jk 


rn jk) + I 1 


l )k 


a 


jk)*k( n - m jk ) 


+ «?***(» - - !) + 


rn 


)k 


+ C'jfc { ( 1 - «J)$ fc (n ~ m+ fc ) + (1 - a yt )$*(n - mj fc ) 


+ <x%$k{n - m jk ~ 1) + - rn jk - 

+ Cj k |(1 - <*£)$*(« - m+ fc ) + (1 - ^^(n - rn jk ) 
+ - m Jfc ~ X ) + a jk®k{n - m~ jk - 


+ 


a 


}k 


n-m+ k ) + (1 - a rt )$ fc ( 


+ a jk$k{n ~ m] k - 1) + a jk $ k ( 


jk> 


n - m jk 


n — m 


jk) 


+ D jk \{ 1 - aj fc )$ fc (n - m+ ) + (1 -- a jk )$ k (n - m jk ) 


+ <*jk®k{ 


+ Djk 


+ Q ]k$k( 


n — rn 


jk 


!) + a jk®k{n - m 


jk 


«;*)♦*( 


n — rn 


jk > 


) + (i -a yfc )$ fc (n- m. fc ) 


jk) 


n — rn 


+ 

jk 


!) + a jk®ki n ~ m 


jk 


where AT has been dropped for notational simplicity (i.e., 4>(nAT) = $(n) ) and 

e % = (™% + a %) AT 


(51) 


(52) 


On the other hand, the time derivative of $ may be evaluated by finite-difference as 

d9 ^ T±) = [**(r±) - $k(r± ~ AT)}/ AT ; 

(53) 

= \*k(r±) - 2 <Mr ± - AT) + $ k (r ± - 2AT)]/Ar 2 


dr 

8 2 mr±) 


dr 2 


where r± = t 0 ± R' 0 and AT is an arbitrary small time increment. 
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Note that if the time increment is sufficiently large, the integer delay-time index (for 
instance, mj k ) might be zero so that at the present time step an unknown would appear 
on the right hand side of Equation (51). The co-normalwash on the other hand is known 
(prescribed) at all times — including the current time step. 

In order to express $*.(n) explicitly (in terms of $&(n — mj k ) ), re-define the integral 
coefficients so that 


B){ = (1 - -t + B jk j AT + B, k /AT 2 ) 

B ',i = (1 - «Jt ) {Bjt + B,„IAT + B jk /AT 2 ) 

B)i = Vv - b , k ( i - 2„;, i/A7' - b , t ( 2 - 
Bj k = Bj k " /k — Bj k (l — 2 a /f ] [AT — Bj k {2 — 3a ]k )/ AT 2 
= ~B, k ^JAT + B jk ( 1 - 3a 2 k )/AT 2 
B% = —Bj k tXj k / AT + /q. i I - 3a jk )/AT 2 
B-i = Bj k ocj k / AT 2 
B] k = B lk a-JAT 2 

' (1 - o&JICyt + (6,4 + B,4)/AT + (C-yt + D ik )/AT 2 } if m+ = 0 
0 if mt ^ 0 


(1 “ <*Jk)[Cjk + {C jk + D jk )/AT + (C jk + Djk)/ AT 2 ] if 


if to . ^ 0 


r i 4 = f (! - + (6;* + + (c jk + b jk )/AT 2 ] if m+ * 0 

I 0 if m j k = 0 

r; i f (1 - «7t)[C >fc + (G> + Djk)/ AT + (C> + I>^)/AT 2 ] if m- 0 
|0 if mj k = 0 

c‘l = CjtaJ - (6,* + £>,*)(! - 2 a+ )/AT - (<?,» + D, k )(2 - 3o+ )/A T 2 
C,* = Cyta"* - (6, t + D, t )(l - Za^/AT - (C, t + X> jt )(2 - 3aT t )/Ar 2 
Cyt + = + P, k )(~a* k IAT) + (C, t + £>, t )(l - 3a+ )/AT 2 

= ( 6,4 + p jk )(- a -JAT) + (c> + - 3 «;j/a r 2 

Cy» = (6,4 + b lk ) a j t /At 2 

c ',i = (6,4 + ftt)«7 t /A( 2 

With equations (54) and (55), Equation (51) becomes 


Y A )k$k{n) = bj(n) 
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with 



Ajk E jk 

s-iO — 

{ 'jk S k 

(57) 

N n 4 

( 



b A n ) = 

- 

■* + 1 ) + B )k 'M n - m jk - i + 1) + 



k 1 i 1 


Cjt ®k{n ~ rnj k - t + 1) + C) k * k [n - m jk —» + !)} 
where Ejk is defined according to Eq. (15). 

Thus, solving the algebraic system of equations (56) yields the solution <S>, at current 
time step n. Note that the right hand side of Equation (56) is completely determined since 
it includes only $ at preceding time steps. 

4. PANEL INTEGRALS FOR THE SOURCE AND DOUBLET COEFFI- 
CIENTS 

The supersonic unsteady Green’s Function numerical formulation, Equation (45), Sec. 
3, involves essentially two types of integrals. The source integrals (Equations (46a~c)) 
possess an integrable singularity on the Mach forecone, so they may be interpreted in a 
classical sense. The doublet integrals (Equations (47a-c)), which involve a derivative of the 
Green’s Function and hence must be viewed in distribution sense in order to be properly 
interpreted. 

With this in mind, we shall formulate the procedure for calculating the supersonic 
coefficients. This procedure is described in detail in Ref. 5. With the exception of the 
two additional integrals arising in the process of introducing the time-domain numerical 
formulation, (Equations 48a and b) the remaining integral calculating procedure is pre- 
sented in Ref. 5. Therefore, for completeness, it is sufficient to give only a summary of 
this procedure here. 

As mentioned earlier, once the coefficients in Equations (45) through (48) are com- 
puted, the potential at time t may be given in terms of the potential and co-normalwash 
at preceding times through the algebraic relation in Equation (56). Due to the complexity 
of these surface integrals, it is impossible to obtain closed-form analytical expression for 
them. As in Ref. 5, we shall follow the tactic of computing these double (surface) integrals, 
in £ and ??, analytically in the £ direction and numerically in the r/ direction. 

4.1 Some Definite Integrals 

In what follows we shall allow the rj integral to be evaluated numerically and will 
analytically compute the £ integral. 

By using equation (29) of Section 3.2, we deduce that 

|gj x a 2 \ = K-Pi + r lP 3 ) x (P 2 + £^3)1 (58) 

We now make the approximation 

|«i(h) X g 2 (£)| ~ A(rj) + B(r?)£ (59) 
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where 


A(ri) = | a.j (rj) X a 2 (0)| 


and 


Bin) 


|ai(h) x «2 (1) | - |°i (h) X a 2 (—l)| 


(60) 


(61) 


Making use of Equation (59) and recalling the relationship (27), the integrals (46) 
through (48) may be reduced to the consideration of the following ^-integrals: 



+ 1 

im{v) = / rHR'di 

and 

KM = f\ t m Hdt 

for m — 0, 1, 2 

The /3 m (??), Imiv) and K m (r?) are convergent integrals while the ‘p.f.’ in the a m {rj) 
and A m (rj) integrals indicate that these integrals must be interpreted in the sense of the 
Hadamard Finite Part in order to assume a finite value. 


4.2 Some Indefinite Integrals 

In this section we shall explicitly obtain the indefinite integrals 

OmU,») = / r| I 
&»<{.») = / r§ 

J rn'd( 

i.(u) = / r|| 

and 

«(e,»?) = J rdn 


(63) 
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for m — 0, 1, 2 


which occur in the process of evaluating the first order supersonic source and doublet 
coefficients. 

For convenience we write 


where 


We also let 


R'=\\R\\ = („e 4 be 4- r) l/2 


a = a\ o aj 

b = 2 (P 0 + ?7 -'P 2 ) ° ai 

c = (P 0 + ’nP 2 ) o (P 0 + rfP 2 ) 

d = b 2 — 4 ac — —4 (J? x aj o x aj) 


(64) 


(65) 

( 66 ) 


From standard integral tables we obtain 



' -(4a£ + 2b)/dR', 

for d ^ 0 

— l/R'(2a£ + b), 

for d = 0 
a > 0 

^■3/2 

for d, — 0 
a = 0 
c > 0 

0 

otherwise 

— (aR 1 )* 1 — (b/ 2a^d 0 

for a / 0 
d 0 

2R' /b 2 4- (2c/6 2 ) /JZ' 

for a = 0 
d 0 

-(a./?')" 1 - (b/2a)a 0 

for a > 0 
d = 0 

£2 C - 3/2/2 

for a — 0 
b = 0 
c > 0 

0 

otherwise 


(67) 
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A third basic integral we shall need is 


0, 


oU,v) = J 


dt 

R 1 


+ 2a£ + 6) / \/a| / %/a 

for a > 0 
d / 0 

— tan~ ] ^(2a£ + b)/2y/—aR'^j / y/—a 

for a < 0 
d / 0 

ln\2R'y/a + 2 a£ + 6|/-^/a 

for a > 0 
d = 0 

2R'/b 

for a = 0 
d / 0 

e/v^ 

for a = 0 
d = 0, c 


) for otherwise 

In terms of a 0 ,a\ and (3 0 above we are able to express the integrals „ 


«2 


(«,») = / 


( 2 d( 

R> 3 


= < 


0o ~ b& \ ~ c «o)/« 

{2/3)R' 3 -4cR' 2c 2 / R'\/b 3 

£3 c -3/ 2 / s 


for a ^ 0 

for a = 0 

6/ 0 

for a — 0 
b = 0, c > 0 

otherwise 


01 tt> r l) = / 


7o 


(*,»?) = / 



' 2i2' - ,6f 0 6/2a 

for a / 0 


(2/3)R' 3 /b 2 - {2/b 2 )R'c 

for o = 0 

1 

AiV ! 
^ ; 


b / 0 

R' 

ec -»/2/2 

for a = 0 



6 = 0, c > 0 


■ 0 

otherwise 

' (2 a£ • 

■f b)R' /4a -f (4ac — b 2 )(3 0 /8a 

for a / 0 


d / 0 

(2/36) (6^ -f c) 3 / 2 

for a — 0 



d / 0 

\/c£ 


for a = 0 



6 = 0, c > 0 

■ 0 


otherwise 



In terms of j3 o ,0\ and 7 0 we obtain 71 , ft 2 and 72 

[(2/3)1^ - fn 0 \/2a 


7i(e,»?)= / ZR'dt 


\(2/5)R' 5 - (2c / 3) R' 3 ) / b 2 

esfc/2 


for a £ 0 

for a - 0 
ft ^ 0 

for a = 0 
b = 0, c > 0 


(73) 


= / 


72' 


1 


(£/2a)R' - (3ft/4a)/?j - (c/2a)/3 0 
(2/5)(ft* + c) 5 / 2 /ft 3 - (2c/6)& - (c 2 /ft 2 )/? 0 

e 2 /3^- 


otherwise 

for a 0 


(74) 


>- 0 


i2(t,v) = J e#dt 


( (£ - 56/6a)i? ,3 /4a + (5ft 2 - 4ac)7 0 /l6< 
[(2/7) J?' 7 - (4/5)ci? ,s + (2/3)e 2 i? ,3 ]/ft 3 


e 3 v^/s 

0 

Next, consider the basic integral 

e 


for a = 0 
ft ± 0 

for 0 = 0 
5 = 0 
c > 0 

otherwise 
for a/0 

for a = 0 
ft ^ 0 


for a = 0 
ft = 0,c > 0 
otherwise 

(75) 


■< £ '”> = J-b d( = J 


c nd i 

a£ 2 + ft£ + c 


for m ==' 0, 1, 2 


From standard tables we find 


A 0 (£,?7) 


log[(2a£ -f 6 — -v/a) / (2a£ + ft + \/d)}/\/d for a > 0 

(1 > 0 


2 tan 1 \(2a£ + ft)/\/^]/\/— rf 
-2/(2a£ + ft) 
log|ft£ + c|/6 

t! c 


for a > 0 
d < 0 

for a > 0 

d = 0 

for a = 0 
ft ^ 0 

for a = 0 
6 = 0 , 
c ^ 0 


(76) 
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and 


Note also 



(log |a £ 2 + 6 £ + c| - bA 0 )/2a 

for a 7 ^ 0 


Ai(£,»7) = ' 

tlb- 

clog I 6 £ + c|/ 6 2 

for a — 0 
6^0 

(77) 


e/2c 


for a = 0 
b = 0 , c 7 ^ 0 




(£ - 6 A, — cA 0 )/a 

for a 7 ^ 0 


A 2 (£,??) = ' 

1 

C 2 /2b — cA]/b 

1 

for a = 0 
b ^ 0 

(78) 



£73 c 

for a = 0 
6 = 0 




= J d Z = £ 

«! = / £<*£ = £7 2 

*2 = J £ 2 <*£ = £7 3 


(79) 


In sections 4,3 and 4.4 we shall see how the indefinite integrals a m , j3 m 7 m , A m and k m 
given above are utilized in evaluating a m (r}) , /3 m (rj ) , 7 m (??), A m {rj) and /c m (r?) for m — 0,1 
and 2 . 


4.3 EVALUATION OF THE INTEGRALS FOR FULL PANELS 

In the case of a full panel, i.e. one in which {(£,h)| — 1 < £ < +1,-1 < r/ < + 1 } lies 
entirely within the open Mach forecone = {J?|X — X* < 0 and R o R > 0} the evaluation 
of the a m ,/3 rn and 7 m follows easily. In this situation these integrals are convergent and 
the Hadamard Finite part is not needed. 

The indefinite integrals A m and k m given above are to be utilized in eval- 
uating 7 m( 7 ), A m (??) and K m (rj) for m = 0 , 1 , 2 . 

«m(h) = « m (l,»l) - am(-M) 

Pmix) = Pmihri) - fimi-hn) 

Imiv) = 7m(F V) ~ 7m(~F V) (80) 

A m{V ) — A m ( 1 , 7? ) A m ( ~ F V) 

«m(h) = &m( F *?) “ »?) 

for m — 0 , 1, 2 

In this situation the oc m (r]) , /? m (r?) , 7 m (??) , A m (p) and K m (r/) are analytic functions of 
7 for —1 < 77 < -fl, and the numerical computation of the definite integrals involving 
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these functions and appearing in Section 4 may be carried out by a standard numerical 
integration such as Gaussian quadrature. 

4.4 Evaluation of the Integrals for Partial Panels 

In the situation where a panel lies partially within the Mach forecone the evaluation 
of a m {ri), /3 rn (ri), ') m ( r l), A m (r/) and /c m (/?) takes a bit more doing. We note that integrals 
a m , rri = 0, 1, 2 are singular on the Mach cone — {R\X — X* <0, R o R = ()}. Thus, for 
such that R(£ t ,r)) lies on the Mach cone, we must evaluate in accordance with 

the Hadamard Finite Part. We obtain 

p.f. = 0, m = 0, 1 (81) 

p.f. A m (£*,7?)=0, m- 0,1,2 

On the contrary, the integrals f3 m , 7 m and m = 0, 1, 2, are not singular at (£*, rj) so in 
order to calculate the value of these integrals at that point it is enough to plug that point 
into the expressions for these integrals. There are a few provisos however. A problem 
will occur in calculating /?„(£*, 7) from expression (69) in the case a < 0 and d p 0, since 
!?'(£*, 7) = 0. However the identity 


Po{C,n) 


1 tan-h 2 ^ 


\2 y/^dR 1 


sin 


— a 


2 a£ + b 

\fd 


(82) 


for a < 0, d p 0 together with the fact, that 2 a£* -f b = ±d and d > 0 for with !?(£*, 7) 
on the Mach cone shows that 


/?o(£*,h) 


7T 




sgn (2a£* + b) 


for a < 0, dp 0 

In addition, we show, in Appendix A, that 

p.f. d 2 (£*,7) = { 2( Z a) 3 ^ Sgn ( 2a t* + 


(83) 


for a < 0, d = 0 
otherwise 


(84) 


Since the expressions /?i,7o,7i,/?2 and 72 are all given in terms of a 0 ,ai,/2 o and R\ 
there is no difficulty in evaluating these functions at £* with #(£*,7) on the Mach cone 
via use of Equations (70) — (84) where applicable. 

We are now ready to investigate our integrals for a fixed 7 with — 1 < r) < 1. We 
look at the intersection of the interval — 1 < £ < +1 with the Mach forecone {R\R o R > 
0, X — X* < 0}. Four cases may occur: 


23 



(i) The intersection is empty. 

(ii) The intersection is a closed interval [£*, £ tt ] with £i < £ u . 

(See Figure 3.) 

(iii) The intersection is a single point £ 0 with |£ 0 | = 1 but R(£ 0 ,p) °oi ^ 0. 

(See Figure 4.) 

(iv) The intersection is a single point £ 0 with 

(a) — 1 < £ 0 < 1 and R(£„,p) ° a.\ = 0 (See Figure 5a.) 

(b) |£„| = 1 and i?(£ 0 ,r/) o a, = 0 (See Figure 5b.) 

A point (£ 0 , p 0 ) where X — X* < 0, R o R = 0 and R oa\ = 0 is called a critical point 
if — 1 < £ 0 < 1, — 1 < p 0 < 1, i.e., {£ 0 iPo) is a critical point in Case (iv) above. 

Case (i) 

In this case we define 


^m(h) — 0, — 0) 7>n(h) — 0 

&m{v) = 0 and K m{p) = 0 


Case (ii) 

In this case we define 

<*m(h) = d w (£ u , p) - a m {£ ( ,p) 

0m{p) = Pm{£ui P) ~ /^m(Ct) P) 
'lm{p) — '7m( 6u> P) 7 m{£l>P) 
^m{p) X m (£ n ,p) — A m (^,iy) 
K m (??) = k m (£ n ,p) - «m(&,»?) 


(85) 

for m = 0, 1, 2 


(86) 


for m = 0, 1, 2 


where a m , /3 m , 7 m , A m and are given as in Section 5. We remark however that if 
either £ u or £i or both lie on the Mach cone the evaluation of the a m , /3 m , and A m for 
siich £* must follow (81), (83) and (84) of this section where applicable. 


Case (iii) 

In this situation, we have a limiting situation where either £f — > -f 1 or £„ — » — 1 . The 
functions a m (£, p), p) and p m {£, p) are continuous at such a point and thus we find 

that 

<*m(p) = 0 
Pm{p)~ 0 
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1 m(h) = 0 


A m (h) = 0 

and 

«m(h) = 0 

Case (iv) - (£ 0 ,r/ 0 ) is a Critical Point 


(87) 


for m = 0, 1, 2 


Assume at first that -1 < r/ 0 < +1. Then the equation R o R — a£ 2 + b£ + c = 0 
possesses a double root at £„ = — 6/2a with — 1 < £ 0 < 1, provided a < 0. At this point, 
d — 0. (Note: If a — d\ od\ >0, we cannot have a point on the Mach cone with Ro ~d\ = 0 
unless |J?| = 0 identically.) Now in accordance with Eq. (A. 19) of Appendix A, we find 
that in this case (a < 0, d = 0) 

PoM = -£= 

\J—a 

/Mho) = - (88) 

/Mho) = — F=6 o 

\J-a 

while h m (ho) = 0 and K m (ri 0 ) = 0 for m = 0, 1,2. 

The ct m (? 7 ), m = 0, 1,2 behave in a more complicated manner. These expressions may 
be given in the form 

«rn(h) = am g (h) + - ho) for m = 0, 1,2 

where: 

«c Gg (ho) = 0 

«f g (ho) = 0 (89) 

and 

a re S (r) ) £ 

“ 2 Wo)~ ( _ a)3/2 

and the special distributional contribution to 0 t o , a .\ , and 0:2 is given by 6 (17 — rj 0 ) where: 


SP 

a rn 


_-* ar 

| R • a. 1 x a 2 | 


(90) 


for m = 0,1,2 


In the cases where — 1 < rj 0 < 1 and |£ 0 j = 1 the only change with the above is that 
the special contributions «o P , oq P , and aijp are divided in half. 


Thais 


o Sp 

ULrn 


7T 


a n 


2 \R • a ] x a 2 | 


(91) 
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We have not yet considered the situation where |^ 0 | = 1. Here if |£ 0 | < 1 we set 

0 if rj 0 R o a 2 > 0 

— io 1 ™ j\R • 5j x 0 , 2 1 if r) a R o 02 < 0 


sp 


and if both \rj 0 \ = 1 and |£ 0 | = 1 we set 

0 


sp 


if p 0 R o a 2 > 0 


1 io > ' n !\R • a \ x a 2 \ if rj 0 R ° a -2 < 0 


(92) 

for m = 0, 1, 2 

(93) 
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APPENDIX A — CRITICAL POINT AND SPECIAL SINGULARITY 


In this appendix we shall study singular integrals of the form 

P- f - / S ^R<^ = Fs ^ U- 1 ) 

where S(£,r]) is a polynomial in £ of degree < 1. 

We focus on the situation where a = o a, < 0 in the neighborhood of a point 
(£ 0 , Vo)-, — 1 < £ 0 < +1, — 1 < r) 0 < +1 where R o R = 0, R o fij = 0 and X — A* < 0. 

We shall find that this integral (A.l) exists only in distribution sense as a function of 
rj and in fact takes the ‘value’ 


F,M = 


{£()) Vo) 

R ■ 0,\ X ft'2 1 


V -■ r In 


{A. 2) 


Proof: We study in detail the prototype situation of Fig. 6 where R o a 2 > 0 at (£ 0 ,r? 0 ) 
and the Mach cone intersects rj = +1 at two £ values both with |£| < 1. The situation 
with R o g 2 < 0 at (£ 0 , ?y 0 ) with the cone inverted may be handled by analogy. 


We note that for rj i < r) 0 the cone does not interact the line rj = rji so that clearly for such 
Vi, FsM = 0. 

Further for ifo > rj 0 the integral (A.l) has a 3/2-order singularity at two distinct £ 
values on the line v ~ V\- By definition of the Hadamard Finite Part it follows that 
F s {rj ]) = 0 here again. 

In order to demonstrate (A. 2) we therefore need only show that 


J i I R • a .] x a 2 1 


<1 - rio 


In fact if we denote by 


W. 


(v) = P-f. f I 


5(£,P) 

R ri . 


d£ dr] 


n ° allowable ( 


{A- 3) 


(A.4) 


it is clecir from previous remarks about F s {r]) that 


1: 

2 : 


W s {rj) will be independent of rj for r] > r) 0 , and that 



• +1 

( + !) = I F s (rj) dr] 


l 


for V > Vo 
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so that it suffices for us to prove that for any rj > rj 0 

w s [n) = r £ ({ “’ 

» • ». ' «ii ; t 

Interchanging order of integration in (A. 4) we find that for any rj] > rj 0 

r/l 

w»{v i) - p.f. J f —pjP d i 

allowable { r,a 

We now denote by V(i,rj]) the integral 

Vtf,ri i)= p.f. J dr) 

Vo 


We have previously seen that 


jv{i,r))d£ = W s (??) = W s 


RoR>0 
X -X, <0 

is independent of r/ for rj > r) 0 . 

Let us now evaluate the integral (A. 7) explicitly. 
We note that for any fixed rj > rj 0 we may write 


R' - sf- a XiX i)(i - (t) 


(A.S) 


(A. a) 


(A.7) 
for rj j > fj 0 


in the situation of Fig. 6. Here it(rj) and i u (rj) are roots of \R\ = 0 with < £„(??). 

If we now write 


U(t,rj) 

l(«. - 0(f - 6)]’/ J 


we find that 


W t = W, 




(A.8) 


' £ " U((,ri)di 

h, [({.-Otr-M ! 1 / 2 

Next making the change of variable £ — (£ u — f l()p/2 4- (£ u + it) {2 in (A. 8) we obtain 


W s = 


M = J 


-1 


u\{ iu - it)pl 2 + {in + it) /2 , v )dp 

(1 -p2)l/2 


(A. 9) 
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From (A. 7) we have seen that for rj > t) 0 the integral (A. 9) is independent of r?. Now 
letting r] j rj 0 we of course have £ u (^) j £ 0 and (t(r)) j £ 0 and then obtain that 


W c 


l 

u{L,Vo) / (1 _* 2) 1 /2 = *U(to,rio) 


(A. 10) 


Next we compute U(£ 0 , r/ 0 ) in terms of S(( 0 , rj 0 ). 

We note that 

so setting £ = £„ and using L’Hopital’s Rule as 77 j r\ 0 we obtain after simplification 

S[Zo,rio) 


_| ai | (Soaj) 


and then from (A. 10) we have 


W. 


StioiVo)* 


(A. 11) 


-|«i|(i? 0 a 2 ) 

Now in the situation as pictured in Fig. 6 R o a 2 > 0. In the situation where the Mach 
cone is inverted, R o a 2 < 0 at (£ 0 , rj 0 ), (see Fig. 7) but r )2 < rj 0 in (A. 6) causing a double 
sign reversal so that we obtain in either case 


3 .f. / 


s{£,v) 

R'Z 


dti 


S(io,r)o) 


1«l | {R 0 a 2 ) 


i=(o 
r} = r)o 


~ Vo) 


(A. 12) 


To complete our proof we note from Appendix B that at a point (£ 0 .r/ 0 ) with R o R, — 0 
and R o ci] = 0 we have that 


| R ■ <i] x a 2 1 = \o.\\\R o a 2 | 


so that finally from (A. 12) 


Av) = P-f- / 






( C 05 »7o) 


|i2 • aj x o 2 | 


- io 
rt-no 


Tr«5(?7 - »7 0 ) 


(A. 13) 


(A. 14) 


which completes the proof of (A. 2). 


We are now prepared to relate the foregoing to the evaluation of the a m (r]) for m — 
0,1,2 at a critical point (£, 0 iVo)- For simplicity we assume — 1 < £ 0 < 1 and — 1 < r] 0 < 1 
in our discussion. 

We recall that the a m {v) are defined as - 

£ m 


OV 


(v) = p.f. I Ri: , for m = 0,1,2 
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In the cases m = 0 or m = 1, £ m will of course be a polynomial of degree < 1 so that 
from item (A. 2) we can immediately say that with S(£ 0 ,r) 0 ) = ££* 


«m(b) 


&‘Tr 


\ R - a l X «2| 14 = 4° 


Mb - bo) 


In the notation of Section 7 we set a^f^(bo) = 0 and 


a Sp - — 


R • a\ x 0 , 2 ] 


t = to 
ti-rto 


for m = 0 or 1 


for rn = 0 or 1 


(A. 15) 


In the case m = 2we are dealing with 


«2 (b) = P-f. / 


This integral contains both a regular and a singular component. To isolate them let us 
note that if R' = (a£ 2 + 6£ + c) 1 / 2 then £ 2 = R' 2 /a + S(£, rj) where S(f, r?) is a polynomial 
in £ of degree < 1. Note that £ 2 = S(£ 0 ,r) 0 ). 


Thus we may write 




Now the second integral on the right of (A. 16) will equal 


I ^*®1 x s aM 4 = 5 o 


Mb - bo) 


by (A. 2) but since £| = S(£ 0 ,rj 0 ) this integral must equal 

_ -?r£o 

I® -Si x °r 


and we set 

as we have stated in Eq. (83). 
To proceed we must have 


sp _ 


«2 = — 


~7t£o 


|i? • a] x a 2 | A = <j 


a 


reg 


^ a i R' d ^ ■ 


(A. 16) 


(A. 17) 


so 



But this integral equals 


in our notation of Section 7. 
To show 


it suffices to show that 


,Po( r lo) 




0o{*lo) 


(— a ) 3 / 2 


We proceed with the latter. In fact we compute /? m (p) for m — 0, 1,2 at once. Suppose 
the Mach cone cuts our panel as in Fig. 6 and we consider rj > r) 0 . Then with the notation 
R 1 - u _ £)(£ - we see that 


Cti 

Pm{v) - ~F= f 

\ — a .1 


a i y/Uu-m-tt) 


(A. 18) 


for rn — 0, 1, 2 


The change of variable £ = (£ u + £*)/2 + (£„ — £t)p / 2 transforms A. 18 into 


MV) 


i / [(eu + 6)/2 + (eu-6)p/2] f 




dp for m = 0, 1, 2 


V a 4 V 1 “ 

Now letting rj [ rj 0 we have that £ u j £ 0 and £^ f £ 0 so that in the limit 


/ v Z™ [ dp tt£” 

= - 7 = / - 7 = = “ 7 = 

V-«' •' x /l - p 2 V - 


(A.19) 


for m = 0, 1, 2 


An almost identical argument shows that 


'YmW = -Tt) 


for m = 0, 1, 2 


transforms into 


Trn(h) = / [(£« + &)/2 + (£„ - 6)p/2] m [(£„ - 6)/2] 2 /l - p 2 dp 

- 1 

for m = 0, 1, 2 

so that as ?? j p 0 , £ u 1 £ 0 , £< t £o we obtain 

lm{Vo) = 0 form = 0,1,2 (A. 20) 

Finally, consideration of A m (rj) at a critical point yields, in a similar way 

A m (i)„) — 0 for m = 0,1,2 (A. 21) 
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APPENDIX B.— A LEMMA CONCERNING SUPERDOT PRODUCT 

In this brief appendix we prove an elementary lemma concerning the superdot product. 
This lemma comes into play in proving formula (A. 2) of Appendix A. 

LEMMA: Let a, b and c be three vectors in R 3 . Assume that (i) a,oa = Q. (ii) a o b — 0 
and (iii) b o b < 0. Then 

\a • b X c\ — \J — b ob |a o c|. (R.l) 

We refer the reader to Ref. 3 for definition and properties of the superdot product. 


PROOF : Without loss of generality we may assume coordinates have been rotated so that 
a z — 0. We may assume |a| ^ 0. 

Then a takes the form 

a = a x i + a y j with a\ — a} y ^ 0 from (i). 

Let us write 

b = b x i + b y j + b z k 

then 0 = a x b x — a y b y from (ii) or equivalently 

d x b x — ^y b y (B.2) 

We now proceed by cases. 


Case 1: If a x = a y ^ 0 then from (B.2) it follows that b x — b y . Then 

/ 0, x (ly 0 \ 

a ■ h x c = Det b x b y b z 

\ G Cy J 

Q-x{byCz b z Cy) — a y {b x c z — b z c x ) 

Thus a - b Xc = (a o c)b z and \a • b X c\ — \a o c\\J —b o b in this case. 


Case 2: Here a x = —a y ^ 0. 

It follows from (B.2) that b x — —b y 
Then again 

a -b x c = Det 



b, 

c. 
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— ttxipy^z bz c y) a yipx c z b z c x ) 

— (fljCj — 0,yCy ) b z 

= (flo c){-b z ) 

and so ^ 

\a • b x c| = | a o c\\J—b o b once again. 

APPENDIX C- DERIVATION OF THE SUPERSONIC 

OSCILLATORY P.D.E. 

We begin with the linearized potential flow equation: 




0 


Passing to Prandtl-Glauert coordinates after the introduction of scaled 
indicated in equation (3) we proceed as follows: 


1 0 2 4> _ 2 Poo d 2 <\> _ U%o <) 2 4> __ 0 
al o Ot 2 alo dtdx Ox 2 

2 0 2 <f> d 2 4> t 0 V 1 0 2 (\> 2Moo 0 2 (j> 

P Ox 2 + <V aX <9* 2 Ooo df “ 


2 d 2 $ d 2 $ d 2 $ 

~~ * dlJlX] 2 + d(lY ) 2 + 5(FZ) 2 
1 3 2 $ 2M*, <9 2 $ 

0(lT/a^3) 2 ~ d(lT / aoo fi)d(filX) ~ 

which results after simplification in 

d 2 $ d 2 $ 0 2 $ n9 0 2 $ w d 2 # 

OX 2 OY 2 OZ 2 P OT 2 °°OXOT ~ 


Then equation (C.2) can be written as 


0 d 2 $ (9 2 # 

V o V# + 0 2 . )T2 + 2M :x .. )X . )T = 0 


which is Eq. (4) of the text. 


(C.l) 
variables as 


(C.2) 


(C.3) 
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Fig. 1 Geometry of the hyperboloidal element 



Fig. 2 Surface geometry 
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